Confinement-enhanced spin relaxation for electron ensembles in large quantum dots 
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We present a numerical study of spin relaxation in a semiclassical electron ensemble in a large 
ballistic quantum dot. The dot is defined in a GaAs/AlGaAs heterojunction system with a two- 
dimensional electron gas, and relaxation occurs due to Dresselhaus and Rashba spin orbit interaction. 
We find that confinement in a micronscale dot can result in strongly enhanced relaxation with 
respect to a fi:ee two-dimensional electron ensemble, contrary to the established result that strong 
confinement or frequent momentum scattering reduces relaxation. This effect occurs when the size 
of the system is on the order of the spin precession length, but smaller than the mean free path. 

PACS numbers: 72.25.Rb, 73.63.Kv, 73.23.-b 



Due to spin-orbit interaction (SOI), the state of elec- 
tron spins is influenced by electron transport in electronic 
devices. This has been recognized as a source for de- 
phasing and relaxation for spins^i^, as well as a means 
for controlled spin manipulation^ in research that aims 
at developing spintronic devices^i^i^. In this article we 
present a numerical study of spin relaxation in an elec- 
tron ensemble that is scattering inside a micronscale de- 
vice structure. We are interested in the case where de- 
vices are made of clean semiconductor heterostructures 
and studied at low temperatures. For our studies we 
assume realistic material parameters for a system with a 
two-dimensional electron gas (2DEG) at a GaAs/AlGaAs 
heterojunction. For a free 2DEG in these materials, a 
well established result is that the average spin orienta- 
tion of an ensemble decays due to precession in spin-orbit 
fields. For moderate electron mobilities, this so-called 
D'yakonov-Perel' (DP) mechanism for spin relaxation^ 
has the property that the spin relaxation time Ti in- 
creases when the mobility (and thereby the time scale 
for elastic momentum scattering) decreases, as Ti oc t'^^ . 

A similar trend is observed when the degree of electron 
confinement in a device structure is increased. Stronger 
confinement gives more frequent scattering on the bound- 
aries of the system, and hence reduces the relaxation. 
This has been recognized in the increase of Ti on the tran- 
sition from 2D to ID systems^. In the limit quantum con- 
finement in extremely small devices (much smaller sys- 
tems than we consider for the present study), relaxation 
and dephasing due to SOI is then strongly reduced, and 
other relaxation mechanisms can become dominant. This 
applies for example to spin dephasing in few-electron 
quantum dots, which can be mainly due to interaction 
with nuclear spins^'i^'^. Although more frequent scat- 
tering due to stronger confinement thus seems similar 
to reducing the mobility in bulk materials, the results 
that we present here show that frequent scattering due 
to confinement can also result in the opposite, namely 
confinement-enhanced relaxation. 

The key result of the present study is well presented 
by the traces for spin relaxation time Ti as a function of 
the size L of a square quantum dot in Fig. [T^. The trace 
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FIG. 1: (a) Spin ensemble relaxation in a quantum dot system 
of size L. The relaxation time T\ is calculated as a function 
the size L for a square system. For mobility ^ = 100 m^/Vs 
we plot Ti for zero external magnetic field (filled black sym- 
bols) and for B^^t = 10 T (open symbols). The gray symbols 
show Ti for zero external mag netic field and = 1000 mVVs. 
(b) Relaxation time as a function of L for a system with only 
Rashba SOI. Calculations for specular and non-specular re- 
flections give qualitatively the same results, but the magni- 
tude of the resonant structure in the traces of Ti as a function 
of L is larger in the case of specular reflections. 



for the case that an external magnetic field B^xt = T 
and 2DEG mobility /i = 100 m^/Vs, shows that Ti is 
constant for L > 10 /xm. Here L is so large that the 
electron ensemble behaves as in a free 2DEG. When de- 
creasing L below 10 /im, Ti increases because spin re- 
laxation in suppressed by more frequent scattering on 
the edge of the system. The trace for Bext = 10 T and 
2DEG mobility fj, = 100 m^/Vs (typical parameters for 
research on spin effects in micronscale quantum dots^ 
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FIG. 2: (a) Schematic representation of the direction and 
magnitude of the SOI fields. Since we assume all electrons 
move with the same magnitude for fc- vector kp , we can repre- 
sent the motion of electrons in all directions as a circle in the 
{kx,ky)-p\a.ne. The arrows that are sketched at certain points 
on this Fermi circle, represent the strength and direction of 
the SO field Bgo for that fc-vector. (b) Magnitude and direc- 
tion of the total effective magnetic field when an external field 
B^xt = 1.5 T is applied along the [100]-direction. (c) Idem 
for an external field Bext ~ 10 T. In this plot the length of 
all arrows has been scaled down by a factor 5 as compared to 
(a) and (b). 



and wires^^), however, shows radically different behav- 
ior. Now Ti first slowly decreases when decreasing L from 
the 2D regime (very large L), and shows a pronomiced 
dip for L « 1 /xm. Confinement now strongly enhances 
relaxation, instead of the more familiar result that con- 
finement reduces relaxation. Moreover, in the range with 
1 /im < L < 10 /im, Ti has a highly structured depen- 
dence on L. Only when decreasing L below L Pi 1 /im, 
Ti shows again a strong and monotonic increase as for 
the confinement-suppressed relaxation in -Bext = T. 

We obtain these results with a numerical Monte Carlo 
approach. This has the advantage that we can study 
precessional relaxation for realistic conditions, where the 
total magnetic field is the sum of several spin-orbit contri- 
butions and an external field. In reality, samples typically 
have both Rashba and Dresselhaus SOI that are compa- 
rable in magnitude^'*, and in experiments one often needs 
to apply strong external magnetic fields for realizing spin 
transport in non-magnetic semiconductor a^^'^'^i"'^^ . Ear- 
lier studies of these relaxation phenomena were often 
restricted to more tractable cases, as for example with 
only Rashba SOI^^ (no Dresselhaus SOI), and no ex- 
ternal fields. Below, we will summarize our numerical 
method, and then focus on studying the dependence of 
Ti on the degree of confinement in micronscale quantum 
dots. We will also show that in regimes and for param- 
eters that were studied before, our simulations give the 
conventional results. 

We use a description where spin-orbit (SO) coupfing 
acts as a /c-vector dependent effective magnetic field on 
the electron spins. In a 2DEG it is dominated by two 
sourcesl"^. The first arises due to the inversion asym- 



metry in the potential profile of the heterostructure and 
results in an effective Rashba magnetic field B /fii. The 
second effect arises due to the lack of inversion symmetry 
in the GaAs crystal lattice, which is of the zinc-blende 
type, and yields for a 2DEG the linear and cubic Dres- 
selhaus fields Bci and B^iaiJ^. The effective SO field 
Bso in a 2DEG can thus be described as the vector sum 
of these three components, given by^ii^ 

Bfl = Cfl(xfcj^ - yfc^), (1) 
Bdi = CDi{-i^k^+yky), (2) 
Bds = CD3{^ka;ky -ykykl), (3) 

where Cr, Cdi and Cd3 are the coupling parameters, 
X is the unit vector in the [100]-direction, and y in the 
[010]-direction. Our results are calculated using the SO 
parameters that were reported by Miller et al}"^, Cr ~ 
-1.96 • 10-8 Tm, Cdi = -1-57 ■ IQ-** Tm, and Cos = 

-1.18 • 10^24 

The total SO field is then anisotropic in momentum 
space as is shown in Fig. [2^. Each time an electron 
scatters and its direction of motion changes it will pre- 
cess around a different axis set by B^o- For our set of 
SO parameters these fields B50 lie more or less parallel 
to the [110]-direction for almost all fc-directions. Fig- 
ures [2Jd,c show the effect of adding an external magnetic 
field Bext, which is independent of momentum direction. 
We consider here the situation that 'Bext \\ x. The to- 
tal effective magnetic field Bjot is then the vector sum 
of the SO fields and the external magnetic field. When 
the magnitude B^xt of Bext is comparable to that of the 
SO fields, as shown in Fig. for Bext — 1-5 T, the to- 
tal effective magnetic fields are no longer mainly parallel 
to the [110]-direction and there is larger spread in the 
directions of the precession axes Btot- Figure [2t shows 
that for very large external magnetic fields, shown here 
for Bext = 10 T, the total effective magnetic fields Btot 
align with Bext along the [100]-direction. This again re- 
duces the spread in the direction of the precession axes 
Btot- 

In our numerical approach we use a classical descrip- 
tion of the electron motion, and a quantum mechanical 
description of the dynamics of the electron spin. We thus 
assume that electrons have at all times a well-defined 
fc-vector, and electrons move along classical trajectories 
with specular scattering on the boundaries of the system, 
and scattering in a random direction on static potential 
fiuctuations due to impurities. Electrons never escape 
from the system. Although in micronscale quantum dots 
the motion of electrons is confined, the mean level spac- 
ing is much smaller than temperature. Am <S fcsT', 
and this allows for this semiclassical description^^. Fur- 
ther, we consider the case that all the electrons that carry 
the spin orientation are near the Fermi level. Thus, we 
assume that all electron always move with the Fermi ve- 
locity (fc-vectors with magnitude kp), independent of the 
momentum direction. This is a valid approximation for 
ksT, AEz.so ^ Bp (with respect to the bottom of the 
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conduction band), where lS.Ez,so is the Zeeman sphtting 
due to the SO field alone. Obviously, the validity of our 
approach breaks down in the limit of very small dots, 
where electrons are highly localized due to quantum con- 
finement. In practice, this occurs for quantum dots with 
size L below ~ 400 nm, but our results for this regime 
always show a very strong suppression of precessional 
relaxation, which is the semiclassical equivalent for sup- 
pressed relaxation for quantum confined electronsiSiii. 

Our simulation then works as follows. It starts at 
t = Q with each electron at a random position in a square 
shaped dot of size L, and with a fc- vector in a random 
direction. We always consider the case that at t = the 
spin state is prepared in the positive x-direction. For 
each electron, we follow its state in time, and its mo- 
mentum direction will change at each scattering event. 
During each ballistic trajectory between scatter events, 
the electron has a well-defined fc-vector, and we calculate 
the effective spin-orbit field during this trajectory with 
Eqs. [T][31 After each scattering event the electron will 
thus precess around a new effective magnetic field, and 
we follow the quantum mechanical spin evolution in the 
total effective magnetic field during each ballistic trajec- 
tory. We thus find the spin state of each electron as a 
function of time. 

The spin evolution is calculated for an ensemble of 
(at least) 10^ electrons. This mimics the averaging over 
many electrons in an electron transport experiment, since 
large quantum dots behave in practice as a chaotic bal- 
listic cavit y^^'^" . In our model the magnitude of the av- 
erage spin orientation for the ensemble decays to zero 
because each electron has its own scattering trajectory, 
such that the relative difference between the precessional 
dynamics of individual electrons increases in time. We 
will always refer to this as spin relaxation for the ensem- 
ble (with decay time Ti), rather than dephasing, because 
we concentrate on the loss of average spin orientation 
in the direction of the external magnetic field that we 
apply. Note however, that in our description each in- 
dividual electron always keeps precessing coherently in 
its particular spin-orbit field. Consequently, the under- 
lying mechanism is equivalent to that of spin dephasing 
for an ensemble, and in particular for our simulations 
with Bext = T one could argue that the loss of spin 
orientation should be named dephasing. 

The average spin orientation is calculated for the whole 
ensemble as a function of time, independent of the posi- 
tion of the individual electrons. In this study we con- 
centrate on the average spin polarization {Sx) in the 
x-direction. We find in all cases that we consider here 
that the decay time for {Sx) equals that of (S*), where 
{S) — yj {SxY + {Sy)"^ + {SzY, because no significant 
polarization develops in the y or z-direction. The re- 
laxation time T\ is then defined as the time when {S^^ 
is reduced to 1/e of its initial value at i = 0. We used 
electron density 1.0 • 10^^ vcT^ and mobility 100 m^/Vs, 
unless stated otherwise. We neglect inelastic scattering 
mechanisms and electron-electron interactions. 



The momentum direction for an electron changes af- 
ter specular scattering on the boundary of the system. 
This process we will refer to as edge scattering and the 
typical length scale involved here is L, where the area 
of the quantum dot system \s A = 1? . A second effect 
causing a change in momentum direction is scattering 
on static fluctuations in the potential due to impurities. 
Here, a scatter event changes the momentum into a ran- 
dom direction. We incorporate this into our modeling as 
follows. When an electron is moving ballistically through 
the system, the probability that it did not scatter due to 
impurities decreases as e~*/'^= (where Tg is the average 
impurity scatter time) and this probability is reset to 1 
after each impurity scatter event. We thus define the 
mean free path Lm/p = |fF|'''s as the length in between 
scatter events when only considering impurity scattering. 

Another important length scale in our system is the 
so-called precession length {Lpr). This is defined as the 
length of the trajectory where a spin has precessed over 
an angle tt. In our system this length scale is fc-vector de- 
pendent due to the anisotropy of SO fields. Therefore, we 
define {Lpr) as the length of the trajectory for precession 
over an angle tt around the average total effective mag- 
netic field (|Btot|), where we average over all fc-directions 
to account for the anisotropy for the SO contribution to 
the total field. 

Evaluating the relative size of these three length scales, 
L, Lmfp, and {Lpr), for a specific system helps to under- 
stand many properties of the relaxation time. Before 
discussing results for quantum dots, it is instructive to 
discuss two regimes that occur for a free 2DEG {L very 
large). For such systems. Fig. [3] presents traces with Ti 
as a function of Lmfp- We first focus on the case with 
Bext = T, for which {Lpr) = 8 /im. If L^fp < {Lpr), 
the precession angle in between scatter events is small. 
Consequently, the spin state only slowly diffuses away 
from its initial direction in a random walk-like process. 
Here scattering suppresses relaxation. This regime is 
known as the motional narrowing regime^, and the re- 
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FIG. 3: Relaxation time Ti as a function of the mean free path 
Lmfp for a free 2DEG. The top axis shows the corresponding 
value for mobility /i. Calculations for zero external field and 
Be.t = 10 T. 
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laxation time is inversely proportional to the scatter time, 
Ti oc T~^. When Lm/p ^ {Lpr), the spins will coher- 
ently precess over at least one full rotation between scat- 
ter events. Without momentum scattering, the ensem- 
ble shows spin relaxation since each electrons processes 
in a different field, but each electron will maintain its 
component in the direction of its precession axis. Only 
when a scatter event occurs this coherent precession is 
disturbed, such that further relaxation for the ensemble 
can occur. In this regime Ti oc r^. The crossover be- 
tween these two regimes (where Ti shows a minimum) 
occurs for L^fp ~ (Lpr). Switching on an external mag- 
netic field of 10 T for this system, reduces the length scale 
{Lpr) to 1.1 /im. For Ti as a function of Lmfp, this only 
results in a shift of the entire curve, with the minimum 
now occurring at the new value where Lmfp ~ {Lpr). We 
conclude here that our calculations for Fig. [3] reproduce 
the conventional result^. 

We now turn to discussing results for quantum dots, for 
which L can be smaller than {Lpr) and Lmfp- Figure [T^ 
shows simulations of the relaxation time as a function of 
dot size L. We first discuss the result for zero external 
magnetic field (which gives {Lpr) = 8 /im) and mobility 
H = 100 m^/Vs (corresponding to Lmfp — 5 /im). For 
this system, a decrease of L in the regime with L < {Lpr) 
results in higher Ti values. We will refer to this regime as 
the quasi-OD regime. Here frequent scattering suppresses 
precession, as for motional narrowing in the 2D regime. 
Here, we observe that Ti oc a dependence on L that 
was also found for studies on the suppressed relaxation 
in long quasi-lD channels of width L^'-^-. For L > {Lpr) 
we do not observe any change in Ti . This will be denoted 
as the 2D regime, for which L ^ {Lpr), Lmfp- 

For Bext = 10 T (and again /i = 100 m^/Vs) the be- 
havior is dramatically different (also shown in Fig. [T]a,). 
When coming from the 2D regime, there is no longer sim- 
ply an increase in Ti when lowering L towards the regime 
where L <C {Lpr), Lmfp- Instead, there is a regime, here 
for 0.6 /im < L < 10 /xm, where Ti is strongly sup- 
pressed. Moreover, the decrease in Ti when lowering L 
from 10 //TO to 1 /ZTTi shows a structured pattern. For this 
value of the external magnetic field the precession length 
is reduced to {Lpr) — 1.1 /im. The mean free path is 
still Lmfp = 5 /im, such that switching on a strong field 
opens up a regime with {Lpr) ^ L < Lmfp in between 
the 2D and the quasi-OD regimes. Notably, switching on 
a strong magnetic field for a system in the 2D regime 
increases Ti by about one order of magnitude. However, 
switching on a field for a system in the same material 
with L « 1 /im causes Ti to go down more than 2 orders 
of magnitude. Now confinement enhances relaxation. 

For Bext = T and pi ~ 100 m^/Vs we do not see a 
dependence of Ti on L in the regime where L > {Lpr), 
because for that system Lmfp ~ {Lpr)- However, also 
in zero external magnetic field we can open up a regime 
where {Lpr) ^ L < Lmfp by choosing a higher value 
for mobility. This is demonstrated for /i = 1000 m^/Vs 
(which gives Lmfp = 50 /im) in Fig. [1^. Now for 



7 /xm < L < 100 /im the relaxation time decreases when 
decreasing L, and again Ti shows a structured pattern 
(i.e. the structure on this trace here is not noise from 
averaging over a finite ensemble) . 

The reduction in Ti due to stronger confinement is thus 
a general effect and appears whenever {Lpr) ~ L < Lmfp- 
These are typical conditions for micronscale quantum 
dots when large external magnetic fields are applied. No- 
tably, the values that we obtain here for Ti are very close 
to the values that we recently observed in spin accumu- 
lation experiments in micronscale quantum dots^^. We 
found Ti w 300 ps for a quantum dot with L w 1.1 /tm 
and Bext — 8.5 T. This indicates that our simulations 
generate realistic numbers. 

We will now analyze the relaxation mechanism for this 
confinement-enhanced relaxation. To study why there 
is strong dip and structure in the dependence of Ti on 
L, we choose a simplified model system where we only 
consider the Rashba SOI and an external magnetic field 
Bext — 10 T. Both Dresselhaus SOI contributions have 
been set to zero. Using only Rashba has the advantage 
that the magnitude for the SO fields is identical for all 
fc- vectors. The structure in Ti appears more regular here 
(see Fig.lTJj). This proves that the effects that we present 
here do not only occur for very particular SO parameters. 

We have repeated this calculation where we pro- 
grammed non-specular reflections on the walls of the 
quantum dot, such that after hitting a side of the quan- 
tum dot the electron is reflected with random angle back 
into the quantum dot (open symbols in Fig. [TJa). The 
structure on Ti still appears, which confirms that self- 
repeating patterns are not the origin of this effect. We 
observe, however, that the amplitude of the structure on 
Ti is reduced for this setting. This is caused by an in- 
creased variation in the length of trajectories in between 
scatter events for non-specular reflections. We checked 
this by making histograms of the trajectory lengths for 
specular and non-specular reflections (not shown). 

For both traces in Fig. [T|d, we find that the minimums 
in Ti appear at odd multiples of the average precession 
length {Lpr), and local maximums occur at even multi- 
ples of {Lpr)- This means that for systems with a size 
equal to an odd multiple of {Lpr) the electrons scatter (on 
average) after processing (again, on average) over an an- 
gle of exactly tt (mod 2tt) between scatter events, furthest 
away from their original state. This causes fast relax- 
ation. For systems with a size equal to an even multiple 
of {Lpr), the electrons scatter on average after process- 
ing 27r (mod 27r) , so when they are back in their original 
state. Then relaxation is slower as compared to the local 
minimums. However, note that for these local maximums 
there is overall still a reduction of Ti due to confinement 
as compared to free 2DEG: in this regime, more frequent 
scattering on the edge of the system always enhances 
relaxation. The character of the relaxation mechanism 
itself is thus similar to the regime with {Lpr) <C Lmfp 
for 2D systems, where Ti oc r^. The additional feature 
here is the structure on Ti as a function of L, which 
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signals that the overall relaxation mechanism is either 
somewhat resonantly enhanced or suppressed when the 
time-of-flight across the dot matches even or odd multi- 
ples of the spin precession time for an angle tt. Notably, 
the results for free 2DEG in the regime with Ti oc r, 
(Fig. [3]) do not show structure on Ti because there is a 
larger spread in the scatter times Tg. We checked that 
when we program that impurity scattering in a random 
direction always occurs after a fixed time Tg (no spread) , 
we also observe structure on Ti as a function of L^fp for 
a free 2DEG (not shown). 

The most extreme suppression of Ti due to confine- 
ment is in ballistic quantum dots (L <C ^m/p) when 
the size of the system L = (Lpr). Figure [2fc indicates 
that this is a counter-intuitive result when this condition 
is met in strong external fields. The various precession 
axes get more and more aligned when the external field is 
increased to 10 T in x-direction, while the spins are pre- 
pared in this direction. Nevertheless, the lowest Ti value 
that occurs for the various traces in Fig. [T]a, is for a quan- 
tum dot system of L — 1.1 /im in a field of i?ej,t = 10 T. 
Due to the initial spin state in the x-direction, spins 
are initially precessing with relatively small cone angles 
around these effective magnetic fields, and all electrons 
will maintain a large component in the x-direction. This 
results only in a small reduction of (Sx)- Further re- 
laxation for the ensemble only progresses when spins hop 
onto wider precession cone angles, which only occurs at a 
scatter event. Thus, more scattering leads to more rapid 
relaxation, in particular when L is an odd multiple of 
{Lpr), and most rapidly when L = {Lpr). The reason 
that this results in a very fast relaxation mechanism for 
small systems in strong magnetic fields is that the pre- 
cession and scatter times that underlie this mechanism 



are then very short. 

It is interesting to note that a similar conclusion was 
reached with a very different approach in work that stud- 
ied how conductance fluctuations of large quantum dots 
are influenced by spin-orbit effects and strong in-plane 
magnetic fields. Here, it was found that applying a strong 
in-plane magnetic field can enhance the suppression of 
conductance fluctuations by spin-orbit effects^^. Theo- 
retical work on this phenomenon reached the conclusion 
that this effect is strongest in ballistic dots, where the 
precession time is on the order of the time of flight across 
the dot^^-. 

In conclusion, we have shown that the relaxation time 
Ti for a spin population in a certain 2DEG material can 
be strongly decreased when bringing the size of the sys- 
tem from free 2DEG down to micronscale quantum dots. 
The strongest suppression is found for ballistic systems 
of a size L that equals the precession length (Lpr). Here, 
frequent scattering on the edge of the dot rapidly drives 
precession onto wider and wider cone angles, and this ef- 
fect is resonantly enhanced in all systems where L equals 
an odd multiple of (Lpr). We believe our results are very 
useful for comparison to experimental results on this type 
of systems, since we can use realistic device and SO pa- 
rameters. Furthermore, the Ti values that we calculate 
match very well with our recent experimental results on 
micronscale quantum dots^^. 
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